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A mean-field model with a generalized pairing interaction that accounts for neutron-proton pair- 
ing is presented. Both the BCS as well as number-projected solutions of the model are presented. 
For the latter case the Lipkin-Nogami projection technique was extended to encompass the case of 
non-separable proton-neutron systems. The influence of the projection on various pairing phases is 
discussed. In particular, it is shown that number-projection allows for mixing of different pairing 
phases but, simultanously, acts destructively on the proton-neutron correlations. The basic impli- 
cations of proton-neutron pairing correlations on nuclear masses are discussed. It is shown that 
these correlations may provide a natural microscopic explanation of the Wigner energy lacking in 
mean-field models. A possible phase transition from isovector to isoscalar pairing condensate at high 
angular momenta is also discussed. In particular predictions for the dynamical moments of inertia 
for the superdeformed band in ***Ru are given. 
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I. INTRODUCTION 



The strongest effects associated with neutron-proton (np) pairing are expected in N Z nuclei where valence 
protons and neutrons occupy the same shell-model orbits. The basic properties of the np-interaction are known from 
the studies of simple configurations of near closed shell nuclei The isovector (T=l) interaction is dominated 

by the J=0+ channel but the isoscalar (T=0) np-interaction is almost equally attractive in the J=l+ and stretched 
J=(2j)+ channels. The T=0 interaction is on the average stronger than the T=l force. It may, therefore, lead to the 
appearance of a static T=0 np-pairing condensate, particularly in heavier N k, Z nuclei where the large valence space 
allows for the creation of many np-pairs. However, it is not obvious whether these correlations are coherent enough 
to create this new type of collective mode nor what are the main building blocks or specific experimental fingerprints 
of such a condensate. 

Theoretically, np-pairing is a challenging subject. It offers new opportunities to probe specific parts of the effective 
nucleon-nucleon interaction. The generalization of BCS (or HFB) techniques to incorporate and allow for uncon- 
strained interplay of T=0 and T=l pairs on equal footing is by itself non-trivial. Though the first steps to generalize 
the BCS theory as well as the first applications were done already in the sixties |,H,|0||,|l^,|llJ (for review of 
the early efforts see [|l^) only recently the first symmetry unconstrained, self-consistent mean- field calculations have 
been performed [^3|. Extensions beyond mean-field, restoring either rigorously or approximately number symmetry 
and/or isospin symmetry are scarce. 

The renaissance in the interest for np-pairing can be traced back to the fast progress in detection techniques and 
radioactive ion beam (RIB) programs. First experiments with RIB's are soon to come and are targeted on heavy 
proton-rich nuclei in particular on N Z nuclei. They are expected to provide important clues resolving the above 
mentioned, long standing difficulties in understanding np-pairing. The observables to look for are obviously those 
which are expected to be strongly modified by a static np-pair condensate like deuteron-transfer probability Jl^ , p^ , 
/3 and Gamow- Teller decay rates or ground-state and high-spin properties |l^Jl^,|2^,|l|,^ . 

So far no clear, systematic experimental signature of the np-pairing condensate is known. There are, however, some 
indirect indications, for example, in recent spectroscopic data in ggKrag [E3[ and 37Rb37 p4[|. In the ground state 



36^^^ 36 IFJ aii-J 37^ 

of ||Rb37, with T=1,T2=0, the 7— ray energies of the collective 4+ 0+ transitions appear to be similar 

(isobaric analogues) to 3|Kr38, the T=1,T2=1 nucleus in spite of the expected increase in the dynamical moment of 
inertia due to blocking of the like-particle superfluidityE This phenomenon has been interpreted as a manifestation 
of T=l np-pairing collectivity |2^. At higher spins a transition from T=l to T=0 band has also been observed. 
Calculations seem to confirm the T=l np-coUcctivity at low-spins and predict an increasing role of the aligned T=0 



pairs at higher spins, see Refs. [g5 26 2j]. In 3gKr36 a rather unexpected delay of the first crossing frequency has 
been measured ||2^ . It again may have possible links to np-pairing, see discussion in |2^,^,^,^, although more 
conventional explanations involving shape vibrations cannot be ruled out. 

The strongest evidence for the enhancement of np-pairing effects seems to come from binding energies. The well 
known slope discontinuity of the isobaric mass parabola at w Z, see review p9| ] and refs. therein, indicates an 
additional binding energy (Wigner energy) in N k, Z nuclei. The Wigner energy is predominantly due to the T=0 
interaction [pO 31 1. However, the mechanism responsible for the extra binding energy seems to be rather complex 



when expressed in terms of np-pairs of given J,T | |31| . It cannot be solely explained in terms of J=1,T=0 np-pairs, 
at least not for sd or p/ shell nuclei. A connection of the Wigner energy and the T=0 np-pairing condensate was 
suggested in our Letter pO| based on deformed mean-field calculations with a schematic pairing interaction. Mass 
measurements of more heavy N ^ Z nuclei are needed to shed more light on this issue. 

The aim of this paper is to further investigate basic features of np-pairing. The paper supplements the above 
mentioned letter ||20|] explaining in more detail certain technical aspects of our model but also provides new numerical 
and analytical results. The paper is organized as follows: In Sect. |l| we introduce the basic concepts concerning the 
Bogolyubov transformation and the self-consistent symmetries (SCS) used here to simplify the calculations. Details 
concerning the mod el hamiltonian and implications of SCS on the structure and interpretation of the model can be 



found in Section III . Section IV presents the method used to restore approximatively the particle-number symmetry 



which is an extension of the so called Lipkin-Nogami technique for the case of a non-separable proton-neutron system. 
The ideas presented in this section are independent on the kind of two-body interaction used in the calculations. The 
results of numerical calculations, discussion and conclusions are given in Sections ^ and respectively. 



^Throughout the paper, the bold-faced symbols T and Tz=(N-Z)/2 would refer to the total nuclear isospin and its z- 
component, respectively. The T and Tz are reserved to distinguish between various two-body interaction channels. 
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II. THE BOGOLYUBOV TRANSFORMATION AND SELF-CONSISTENT SYMMETRIES 



The starting point of our considerations are the eigenstates of a deformed phenomenological single-particle potential. 
The basis states can be divided into two groups with respect to the signature symmetry {Rx — e^"^" ) quantum number 
r = which are later labeled as ol{cil), respectively. Two different types of nucleonic pairs can therefore be 

formed, namely a.(3 and a.(3 pairs. A generalized BCS (gBCS) theory has to account for a scattering of these two 
types of nucleonic pairs simultaneously. In fact, this work is restricted to pairing in the same and signature reversed 



states i.e. for aa. and aa modes only, see also |11 



The signature symmetry cannot be used as a self-consistent symmetry (SCS) in gBCS calculations. Indeed, in 
such a case the pairing tensor, k — V*U'^, connects only states of opposite signature [|2|. Consequently, the act 
np-pairing cannot be activated, see also ||l^,|l9|. Therefore, to take into account simultaneously ota and aa pairing 
one needs to extend the Bogolyubov transformation. The most general Bogolyubov transformation can be written as: 



a>0 



+ Vakda + C^QfeO- + VakCLa) 



(1) 



where a{a) denote single particle states (including isospin indices) of signature r = —i{+i) respectively, while k denotes 
quasiparticles. As discussed above, by superimposing any SCS one always excludes certain interaction channels in 
the mean-field approximation. Nevertheless, the use of SCS appear many times inherent to the nature of the physical 
problem and of course, make the theory more transparent and easier to handle numerically. Therefore, prior to 
construct a general theory involving the transformation we simplify the problem by superimposing (beyond parity) 
the so called antilinear simplex symmetry, — PTRz, as SCS, see p^ , p^ . One should bear in mind that due to 
the antilinearity of the transformation properties of creation and destruction operators with respect to S'J will 
depend on the phases of the basis states, i.e. no new quantum number can by related directly to this symmetry. 

When superimposing iSj as the SCS it is rather convenient to choose the phase convention in such a way that the 
basis states will have exactly the same transformation properties with respect to both Rx and 5'J and: 



(2) 



Let us divide our quasiparticle states (0) into two families denoted as k and k, respectively: 

al '^{Uakai + Vakaa + Uakal^ + Vakaa) 



(3) 



a>0 



Enforcing S'J symmetry as SCS requires that the quasiparticle operators of eq. (^ have the same transformation 
properties with respect to 5'J as the single particle operators (||) js^. This leads to the following restrictions for the 
coefficients of the Bogolyubov transformation (0): 



f Uak = ] 




ak ak 
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(4) 



The formalism is complex but the density matrix, p = V*V , and pairing tensor k = V*U 
simple structure with real and imaginary blocks decoupled from each other: 



Re{pap) 
Re(p-^) 



iMPap) 
lni{pafj) 



take the relatively 



(5) 



Re(K^^) 
Re{Kai3) 



Im{Kal3) 

Im(K~^) 



(6) 



Furthermore, the complex structure of the single particle potential F and the pairing potential ^ 
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Tq/S = ^ Wa7/3<5P57 Aq/3 = - ^ Wq^^^K^^ (7) 

7(5 7(5 

and consequently the gBCS equations are fully determined by the p and k matrices, respectively. 

In this work we define the two-body np-pairing interaction in terms of an extension of the standard seniority pairing 
interaction. It is separable in the particle-particle channel, Vap^s oc gafisXsi with gap proportional (up to a phase 
factor) to the overlap (qit-I/^t') between the single-particle wave functions .0 As already mentioned it takes essentially 
cxa and cxct. types of pairing, see Sect. |lll| for further detail. To further visualize the physical implications of the iSj 
let us consider the limits of isospin {N = Z case without Coulomb force) and time reversal symmetry (non-rotating 
case). Let us consider ota pairing which in principle consists of both T=l and T=0 components of the Tip-pairing. 
By decomposing the pairing potential into the different isospin components TjT^ one finds: 

i.e. the T=l and T=0 components of np-pairing depend on the combinations of the same elements of the pairing tensor 
but with opposite sign. This is due to the phase relation for the Clebsh-Gordan coefficients which is (anti-)symmetric 
with respect to the interchange of a proton and neutron for T=(0)1, respectively. Time reversal symmetry further 
implies that Q: 



"01/2,0-1/2 '*al/2,a-l/2 '^Q-l/2,al/2 



(9) 



and therefore 



A'^'°^ cxRe(K — -ttt) and A^°'°^ cx Im(K , ,^ — r-rr) (10) 

at^,a-t^ ^ Ql/2,Q-l/2' at^,a-t^ ^ q1/2,q-1/2'' ^ ' 



Consequently, with the pairing tensor of the form of (^, the T=0 component of aa is ruled out through the 5'J 
symmetry. Similar analysis shows that T=l component of the aa pairing also vanishes due to symmetry reasons. 

The latter is well justified due to the Pauli principle. The lack of T=0 aa is a deficiency of our model. However, for 
hu = the aa and aa np-pairing phases are in many applications indistinguishable due to time-reversal symmetry. 
Since our interaction is essentially structureless, based on pair counting we expect our results not to be sensitive to 
this restriction. For hui ^ 0, on the other hand, one wants to associate the T=0, aa (aa) np-pairing with the 
coupling to maximum (minimum) spin, respectively. The retained component is expected to be dominant but only 
at high spins. In order to probe the transition from low to high spin regime one has to explore all possible T=0 
pairs and allow for an unconstrained interplay between the different pairing modes. Note, however, that again due to 
the simplicity of our interaction certain features of the transitional regime can be simulated to some extent with an 
isospin broken hamiltonian. Indeed, the missing T=0 aa component is expected to respond to nuclear rotation in a 
similar way as T=l aa. It was shown explicitely in Ref. p3| for a single j-shell model. 

In conclusion, in our model the aa pairing is equivalent to T=l and aa to T=0 and the isospin notation will be 
used in the following. This simple analysis reveals also the important role the self- consistent symmetries can play in 
theoretical description of rtp-pairing in the mean- field theory, see also pO| . 

III. THE MODEL HAMILTONIAN 

The multipole-multipole expansion offers a rather good approximation to the pairing energy when neutrons and 
protons can be treated separably. In the following we extend this idea to the case of np-pairing by constructing a 
generalized pairing force separable in the particle-particle channel: 

Vpatr = 2^^Vai3jsalalasaj = -^^gapalal o^g*gasa^ (11) 

Q/3 7(5 

where ga/s = {a\G\P) and G is an auxiliary operator generating the specific pairing mode with strength G. The 
antisymmetry of the two-body matrix element Vafj-yS implies that 



^Subscript r in Ut is necessary to distinguish between proton or neutron single-particle states. 
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{a\G\(3) = -{(3\G\a) V a, p. 



(12) 



and therefore, the generators G must be antilinear and antihermitian. We assume here that the correlation energy 
of the nucleonic pair is proportional to the overlap (q;t-|/3t') between single-particle states they occupy (extended 
seniority- type pairing interaction). The G can then be chosen, for example, as: 

rfi /-^aa /^T— 1 rji y^aa /^T— f'l o\ 

for pp{nn) pairing, a.a type of np-pairing and aa of np-pairing. This choice is, however, not unique. In particular, it 
depends on the choice of the relative phases between neutron and proton states. The choice of phase convention induces 
strict transformation rules for the isospin Pauli operators fi, i = x,y,z with respect to time reversal symmetry 
Constructing the generators ( p^ ) we assumed the same phases for proton and neutron states ( \aT,r,T) denotes the 
basis state a of signature r — ±i and isospin r = 1(— 1) for neutrons (protons)): 

f\ar,r^^i,T) = T\ar,r = ±i,T) (14) 
Our phase convention further implies that: 



TT^\ar,r = ±i,T) T\ar,r = ±i, -t) = zf\a^, r = ^i,-T) 



(15) 



leading to Tf^T ^ — t^. Similar considerations for Ty and operators give: 

(16) 

It is straightforward to prove that, with the above phases, the operators ([l3|) satisfy the antilinearity and antihermicity 



requirements formulated in Eq. (12). 

As a main consequence of the separability of the pairing interaction there exists an average gap parameter: 

Using the generators (|l^) we obtain: 

ALt--^-/^^^"^^ -here A^- = G^- ^ (18) 

a-r>0 



{ar\P-r)l^l;' where Aj=i = Ig^^i ^ (a„|/3p) {«^^~ + «^^~} (19) 



Al=o_^=^r{ar\f3^r)IMAl=°) and A^L^iU = -^rK|/3_.)Im(A^=<') 



where A^=o = |Gj=" ^ (a„|/3,) {lm(^~~ ) - Im(«„„^J) 



(20) 



for pp(nn) pairing, T=l, aa type of np-pairing and T=0, aa of np-pairing, respectively. 
The single-particle potential h takes the following form 

hal3 = eaSafj - ^ja(3 + ^ ajS (21) 

where the single-particle energies, Cq,, are calculated using a deformed Woods-Saxon potential ||3^. Nuclear rotation 
is included using the cranking approximation and T originates from the contribution of the pairing interaction to 
the single-particle channel. For pp{nn) pairing we have: 
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r'T=i „ 



(22) 



For T=l, aa np-pairing we obtain: 




Gl;^ J2 {arb-r)(Pr'\S-r') 



Finally, for T=0, cxa mode of Tip-pairing we get: 

1 



pT=0 



ax/3,/ 



K|7-.)(/3r'|<5-r') 




-ilm(p~ 



(23) 



(24) 



IV. THE LIPKIN-NOGAMI METHOD 



The atomic nucleus is a mesoscopic system and as such never undergoes sharp phase transition. The fluctuations 
are always of importance in the transition region. The classical example is the mean-field prediction of a sharp 
superfluid-to-normal phase transition, induced by fast rotation, which is always smeared out in nature. The effect can 
be accounted for theoretically by restoring particle number and this motivated us to include approximate number- 
projection in our model. 

In our letter [pof we have demonstrated that number projection results in a mixing of the T=0 and T=l pairing 
phases already in = Z nuclei. This suggests that the exclusiveness of these modes discussed in the literature is due 
to the mean- field approximation. For example, in the SO (8) model calculations the exact solutions allow for mixing 
of T=0 and T=l phases while not the gBCS |3^Jl^. Recently, it was shown by Goodman pl| , p^ that the mean-field 
calculations with a G-matrix interaction performed in a rich model-space in fact do allow for coexistence of T=0 and 
T=l pairing phases. The number-projection technique used in our calculations is known as the Lipkin-Nogami (LN) 
theory. Below, we will outline the main modifications necessary for applying the LN method to the non-separable 
proton-neutron system. 

The LN theory attempts to construct the state \LN) where both linear and quadratic constraints (tt' G {p, n}): 



{LN\LN)^l, {LN\ANr\LN) 0, (L7V| AiV^ATV^- |L7V) = 



(25) 



are simultaneously fulfilled. Variation over the Lipkin-Nogami state \LN) is equivalent to a restricted HFB-type 
variation for the auxiliary Routhian: 



whereE 



5{LN\H'^\LN) = S{HFB\n'^\HFB) 



(26) 



(27) 



The LN method is not a variational approximation. Only the standard linear constraints for the particle number 

(2) 

are taken into account as Lagrange-type variational constraints. The parameters A^^, are kept constant during the 

lill: 



variational procedure and eventually adjusted self-consistently using three additional subsidiary conditions 

{'H'^iANrANr' - {ANrANr'))) = 0. (28) 

with ANr =Nr-Nr. 



^Obviously, X^p = Ap5 and the symmetric form of the auxiliary Routhian is chosen for convenience. 
*In the following, all averages over \HFB) state will be denoted as ( ) for simplicity. 
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At the point of self-consistency Tijo = O7 ^^'^ the equations ( p8| ) can be rewritten as 

{'H'^\A){A\ANrANr') = (29) 

where |4)(4| denotes the projection onto the 4-quasiparticle space. The LN conditions (^) depend therefore only on 
the two-body residual interaction. These equations form a system of three linear equations: 

J2 AlX'xr"r"' = B^r' (30) 
r"T"' 

wherei 

a;,, = A(.2^) while Xpn = 2X^^2 (31) 

and 

Al'X' = {ANr"ANr"'\A){A\ANrANr') (32) 
Brr' = {Vtu,o-body\'i){MANrANr') (33) 

The x-rr' (or A^^') parameters are therefore equal to: 

(34) 



DetA 



where A"^^ denotes the matrix obtained from matrix A by replacing the rr'-th column of A by the vector B. 

Obviously, the LN theory is technically an analog of the standard HFB theory but for the routhian (p7|). The 
resulting LN equations take therefore the form of HFB equations with the single particle field and pairing field 
renormalized as follows: 

h^^ ^ Kr' + 2\^^lprr' and A:^,^ ^ A^r' - 2\^^},Hrr- (35) 

The contribution to the total energy/routhian due to the LN corrections can be written as: 

5ElN — ~ ^ ^^rT'{PaT,l3T'{5aT,l3T' — PI3t' ,aT) — l^aT,l3T'l^*pr' .ar} (36) 

One should always bear in mind that the LN theory provides a sort of optimal HFB-like wave function being 
different then the \LN) state. Therefore the expectation values of all physical operators should be recalculated using 
the following formula: 

{LN\Q\LN) = (Q) = (Q) - ^ Af^l(g)(A7V,A7V,,) (37) 

tt' 

with \\.^,{Q) dependent on the operator Q in question. These coefficients can be calculated from the following set of 
five equations: 

{QANr) = 0, {Q{ANrANr' - {ANrANr'))) - 0. (38) 



^The asymmetry is induced by the symmetric form of the routhian (p^. Note, however, that the corrections to the F and A 
potentials due to LN terms comes out to be symmetric, see eq. (pq). 
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V. NUMERICAL CALCULATIONS 



An open question in mean-field calculations with rtp-pairing is related to the strength of the interaction. Whereas 
the strength of the pp- and nn-seniority pairing force can be established by a fit to the odd-even mass differences, very 
little is known about the strength of the np-pairing force. Based on isospin invariance arguments, it seems well justified 



to assume that at the N~Z line C'^'^l^) ~ G'^p^. Therefore, we will take in the following G^-^ = {Gf-^ + G^f^)/2 



with G^4p) 



calculated using the average gap method of Ref. ||3^. Concerning the strength of the isoscalar np- 

interaction we will present most of our results as a function of the parameter x^^'^ that scales the strength of T=0 
np-pairing interaction with respect to G^^^ i.e. x'^^'^ — G'^p'^/G'^^^. In some cases, however, we will show examples 
of solutions for the hamiltonian with broken isospin sym metr y. 
The discussion is divided into two parts. 



In section V A we shall concentrate on the solutions for non-rotating 



cases. Subsections VA 1 and V A2 briefly recall major properties of the solutions of both unp rojcc tcd (latcron called 



BCS) and number-projected (lateron called BCSLN) versions of the model. In subsection |V B| , we will discuss a 
procedure allowing to estimate locally, i.e. in relatively narrow range of mass number A, a 'physical' value of the 
scaling parameter a;J^°. In subsection |V C| we shall present the results of cranked, number projected calculations 



demonstrating different aspects of the interplay between mean-field, pairing forces and nuclear rotation. 



A. The frequency zero case 

1. Basic properties of the BCS solutions 

Let us first briefly recall the main properties of the BCS model solutions without number-projection [ pO|jr^ . Let us 
start with the solutions for self-conjugated N — Z nuclei. Disregarding the Coulomb interaction and assuming that 
Gpp"^ — Gnn^ — G-^p^ leads to Aj^f ^ — Aj=^ = Ao. The solutions fall then into three classes depending on the 
value of x'^^^: 

(i) For x"^^^ < 1 {G^p^ < G^p^) the T=l pairing is energetically favored over the T=0 pairing. The pairing energy 
depends only on A'^ = 2Aq -f (A^^^^)^ and no energy is gained by activating the T=l np-pairing. Note that in 
our model all pairing phases can be simultaneously present in contrast to the results presented in Ref. pof . 

{a) The solution at a;-^=° — 1 {Gj-^p^ — Gj^p^) is highly degenerate. The BCS energy depends only on A'^ = 
2Aq -I- (A'^p^)'^ + |A^jj^°p. Also in this limit, no energy is gained due to np-pairing. 

{in) The solution at a;^"° > 1 (G^j^° > G^p^) corresponds to a pure T=0 np-pairing phase. 
For N ^ Z the main features of the BCS solutions can be characterized as follows: 
[w) The T=l, T2=0 np-pairing never forms a collective phase. 

{v) There exists a critical value of the x'^^^ parameter, x'^~^. For x-^^" < x'^~^ only the T=1,|T2| = 1 solution exist 
while for x^^^ > x'^^f the T=1,|T2|=1 and T=0 phases are mixed. 

{vi) The value of x'^f'^ sharply increases with increasing neutron/proton excess jT^j. In consequence there exists a 
critical value |T^''**| above which there is no collective solution for np-pairing. Our calculations show that this 
value is small (jT^*^**! ~ 2), see next subsection. 

Many of the above properties can be understood in terms of the variational approach, inherent to the (generalized) 
BCS-formalism. For the case of equal weight of each pair, the system will always choose only that kind of pairs that 
result in the lowest energy of the system, i.e. result in a sharp phase transition from T=0 to T=l and vice versa. An 
interesting feature of the gBCS solution is the 'redundant' role of the T=l np-pairing in our model when applied to 
even-even nuclei. It does not mix with the other pairing phases and only in N = Z nucleus it coexists with |T2|=1 
pairing. Similar properties were found in Ref. p^ . In the particular case of an N=Z nucleus, the coexistence of T=l 
pairing phases at x^^^=l is very fragile. After switching on the Coulomb force the T=l np-pairing phase essentially 
vanishes. On the other hand, for odd-odd nuclei, one expects the T=l, |T2|=0 pairing to play a role. 

The gBCS solutions for isosymmetric N = Z nuclei have been analyzed analytically in Ref. [|ll[ assuming time- 
reversal invariance and spherical-symmetry in isospace i.e. (Ti)=0 for i = x,y,z. It has been demonstrated that 
coherent, rtp-paired, solutions can be obtained for this case provided that Pap.ap — Pan,an and Kap.ap = —Han.an- The 
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method of Ref. can be generalized to N Z nucleus by extending symmetry condition to (Ti)=0 for i = x,y and 
(T^) = (iV — Z)/2. For N Z nuclei, however, neither Ipap^apl = |pc<n,cm| nor |Kap_Qp| = iKmi.anl- It is then rather 
straightforward to show, using the HFB conditions pn — up* and f? — p ^ kk*, that the gBCS theory does not 
allow for np T=l pairing solutions. For details, see Appendix A. This result relates the strong limitations of possible 
solutions of the generalized BCS theory to the underlying symmetries. In this particular example the limitations are 
due to the time-reversal symmetry and the symmetries ('deformations') in isospace. 

In our model, (Ty)=0, due to the particular form of the density matrix (^. We could verify numerically for 
a number of examples that the BCS solutions for an isospin invariant hamiltonian are realized at (Ta;) = (Tj,)=0 
independent on x^^*^ . This is not surprising since our interaction is symmetric in isospin space. For an isospin broken 
hamiltonian with G^^^ > Gp^^l^^^ and with an overcritical value of G^p^ the solutions for N ^ Z nuclei are triaxial 

in the isospace i.e. with {Tx)j^O, see also discussion in Appendix A. 

One may further extend the gBCS model and include an isospin cranking term, —fiT^. Such a model can be con- 
sidered as the lowest-order approximate isospin-projected gBCS theory. One-dimensional isospin rotations have been 
studied in Ref. in an exactly soluble model including T=l pairing correlations only. It has been demonstrated 
that the isospin-cranking solutions rather poorly approximate the exact solutions but the particle-number projected 
isospin-cranking model offers a very good approximation to the exact solution. Note, that the isospin-cranking model 
may have several interesting analogies to the well-studied cranking approximation for spatial rotations. Due to the 
transformation rules of the isospin operators under time-reversal (Rw which are different than those for ordinary angu- 
lar momentum operators, there may also exist basic differences. Neither the possible analogies nor the differences were 
studied so far. The isospin cranking model in its simplest one-dimensional or more sophisticated three-dimensional 
version, including both T=0 and T=l np-pairing correlations, may provide important clues on the role of np-pairing 
in excited configurations. 



2. Basic properties of the BCSLN solutions 

The BCSLN solution is qualitatively similar to BCS. The major differences can be summarized as follow: 
(i) Both the critical value of a;^^~P as well as the 'physical' value of extracted according to the prescription of 



subsection VB are larger in the BCSLN model than in BCS. For self-conjugated N = Z nuclei, however, the 
ratio of x'^^t / ^o^^ very similar in both BCSLN and BCS models while for 7^ this ratio is smaller in 
number-projected theory. 

(ii) For self-conjugated N = Z nuclei, the solutions of BCSLN allow for mixing of the T=l, jT^j^l pairing phase 
with isoscalar Tip-pairing but not for the collective isovector np-pairing phase. 

A shift of the critical value of x'^^^ towards values that are larger than unity ior N = Z nuclei is due to the asym- 
metric way the Lipkin-Nogami corrections modify the different pairing channels. Whereas the 'normal' constraints 

on proton and neutron number affect the normal and abnormal densities for protons and neutrons, respectively, the 

(2) 

constraint on ANZ involves an additional constrain on pnp (k„p). It appears (numerically) that Xpn is negative 
implying that the Lipkin-Nogami correction due to np-pairing, SE^^ , is positive, i.e. repulsive. Both quantities 
have opposite sign as compared to the analogous quantities for the like-particle channel, see Fig. |l|. It means, that 
the effective gap parameters are weakened in the np- and enhanced in nn- and pp- channels resulting in the above 
mentioned shift, see Eq. (^5|). It is interesting to note, however, that the onset of np-pairing suppresses like-particle 
correlations, increasing both A^^^^^ and the (attractive) energy contribution SE^^J^^"^ in such a way that the total LN 

correction SE^^^^ increases with increasing strength of the isoscalar np-interaction. Note also, that the LN corrections 
and A*^^^ parameters are strongly peaked for = Z nuclei and decrease rapidly with the value of — Z\ (Fig. |l|). 
The A*-^-* parameters can be approximately estimated as: 

(2) 1 d'^E (2) 1 d'^E 

where E{N, Z) is the nuclear binding energy. Using the standard liquid-drop functional form for E{N, Z) one can 
show that the leading order contributions to (^ for A^ « Z nuclei: 

Kin) ^ - + :4 + ^—^N,z (40) 
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x(2) _ las ai ryO-w . , . 

arise from the surface energy, symmetry energy and Wigner term, respectively. Obviously, these estimates are valid for 
complete HFBLN calculations including contributions to the A*^^-* parameters coming from the particle-hole channel. 

(2) (2) 

Nevertheless, the above estimates clearly show that (i) A^^^^ and Xpn have opposite signs and that (ii) an enhancement 

(2) (2) 

in both Apj.^^ and Ap„ is expected for N = Z nuclei due to the singularity of the Wigner energy. These features are 
in nice qualitative agreement with our calculations, see Fig. |l|. Whether or not this behavior is realistic; reflects 
deficiencies of the LN number projection scheme [see discussion in Rcf. |4^] or calls for isospin-projection remains to 
be studied. 



B. Estimate of the isoscalar np-pairing strength 

In our Letter we have demonstrated that the T=0 np-pairing field can yield a microscopic explanation of the 
Wigner energy in even-even nuclei. It naturally accounts for the singularity of nuclear masses in N = Z nuclei. The 
primary mechanism leading to the Wig/-ner cusp can be viewed as a generalization of the well-known blocking effect 
with additional neutrons or protons outside an N = Z core playing a similar role as the odd particle plays in the 
standard blocking phenomenon. This is visualized schematically in Fig. |^. It needs to be stressed that the generalized 
blocking is not a specific property of our schematic model but does apply for a more realistic pairing interactions as 
well, see discussion in [ p7y42| ]. Moreover, while the standard blocking mechanism requires different trial wave functions 
for odd and even systems, the generalized blocking does not. The same trial wave function is used for N = Z and 
N ^ Z nuclei in the case shown in Fig. ||. 

The Wigner cusp can be obtained either by an isospin invariant model with G'^^^ > G"^^^ or isospin broken model 
with G'^p^ > G^^^ppj as shown in ji^. We do not see, however, any particular reason to invoke an isospin-broken 
model, particularly in the light of the discussion in Refs. pO| , pT| . In these works both experimental and theoretical 
arguments are given that the Wigner energy originates essentially from T=0 correlations. ^From these studies it 
is however not clear whether it is due to np-pairing. In fact, in Ref. the np-pairing mechanism was excluded. 
Note, however, that all these studies were performed using the nuclear shell-model. The shell-model Hamiltonian 
is usually written in the particle-particle representation but it can be rewritten in the particle-hole representation 
using the Pandya transformation |^^. Therefore, in the shell-model there is no distinct division into the pairing- and 
single-particle field. This division is inherent to mean-field models only. Indeed, the shell-model definition of pairing 
in terms L = 0, S = 1,T = 1 and L = 0, S' = l,T = Oas derived from the G-matrix interaction in Ref. |^] and used 
in Ref. |Q| to analyze the Wigner energy, is completely inappropriate from the point of mean- field model calculations, 
see also discussion in Ref. pl| , ^ . 

Ref. [0 shows the technique to extract the strength, W^(A), of the Wigner energy [i?vi^(yl) = VF(A)|A^ — Z\] from 
experimental data. Using this method and assuming that the Wigner energy is indeed due to the isoscalar np-pairing 
correlations allows us to determine a 'physical' value of the scaling parameter Xq^^ = G'^p'^/G'^p^ simply by matching 
experimental and calculated values of 14^(^4). Note that this prescription is by no means limited to our schematic 
pairing interaction. Under the assumption that the Wigner energy is indeed due to isoscalar pn-pairing the method 
can be used for any pairing interaction to establish an overall scaling factor between isovector and isoscalar part. 

Let us stress again at this point, that the 'physical' value Xq^'^ deduced in this subsection refers only to the frequency 
zero calculations because of the missing T=0 aa component. For a pairing force, that contains both Cka and aa 
components in the T=0 channel, the same prescription to determine the strength can be used. Note that also for this 
case the same value of G^^° will be obtained. Indeed, due to time reversal invariance, the BCS or BCSLN energy gain 

caused by the T=0 pairing will depend only on th e mo dulus | A-^=°| of the total isoscalar gap, A-^=° = (A^^°, A^g "), 



but not on its direction (see discussion in Section VA). 

Fig. II illustrates calculations for p/-shell A w 48 nuclei. In the calculations we have taken a constant deformation 
of P2 = 0.25 and included 30 deformed neutron and proton states. The realistic values of x^^" are estimated to be 
~1.13 and ~1.30 for BCS and BCSLN models, respectively, as shown in the figure. We performed similar BCSLN 
calculations for A w 76 nuclei assuming /32 — 0.40 and taking 40 deformed neutron and proton states. In this case 
we have taken W^(A) = 47/AMeV ||3l|] as a reference value because experimental masses are not available. The 
calculations yield x^^" ~ 1.25 which would imply that G^^" has a different mass dependence than G^^^. 

In this context one should mention, that the deduced 1/A dependence of VF(A) is essentially due to the fast increase 
of VF(A) for very light sd-shell nuclei. For sd-shell nuclei, however, we expect the T=0 pairing to have a pronounced 
L = 0, S = 1 component. This component becomes strongly quenched due to the spin-orbit interaction when entering 
the p/-shell. The properties of T=0 pairing may change strongly at the border of the pf and sd-shells. The 1/A 
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dependence of may therefore even be an artefact. With the present set of data one cannot extract a rehable 

A dependence of the Wigner energy based only on p/-shell and heavier nuclei. No doubt, that data on masses along 
the N = Z line are urgently required in order to resolve this important question. 

The procedure to determine the strength is very sensitive to even small variations in x'^^^ particularly for BCS 
calculations because we are working in the region of the phase transition. Indeed, changing x'^^^ by ~ 2% (in v4 « 48 

nuclei) affects the Wigner energy by ^ 200 keV, as demonstrated in Fig. |[ Moreover, although a^j^BCSLN ^ ^o|bcs 
in N = Z, A Ki 48 nuclei, the ratio y = Xq^'^ / x'^^^ appears to be very similar in both BCS and BCSLN models. 
However, for ^ we get ?/bcs < J/bcsln as shown in Fig. ^. In consequence iT^ygcgLj^l > IT^IbcsI expected. 
The difference is still rather small and |T^' **| ~ 2 independently on the model. Fig. ^ compares calculated (BCSLN 
model) mass excess AE ~ E{x'^^'^) — E{x'^^'^ — 1) for Cr and Sr isotopes as a function of T^. No pronounced 
differences are seen between Cr and Sr isotopes although there is clear tendency for AE{Tz) to be more spread in 
heavier nuclei, see also Fig. ^. 

The Wigner energy has an additional repulsive component in odd-odd N = Z nuclei. In ||3l| it has been demon- 
strated that the strength of this component d{A)T=Q ~ W^(A) i.e. both components of the Wigner energy have, most 
likely, the same origin. This term can be understood within our pairing theory provided that the odd-spin, T=0 
states in N — Z odd-odd nuclei are not treated on the same footing like even-even nuclei but are interpreted as 
two-quasiparticle {2qp) configurations involving one broken np-pair. Indeed, treating even-even and odd-odd nuclei 
on the same footing within the generalized BCS(LN) approach gives no odd-odd versus even-even mass staggering, 
see Fig. ^. The blocking of a T=0 np quasiparticle, on the other side, would result in an odd-odd versus even-even 
mass (T=0) staggering in N = Z nuclei due to the T=0 pairing energy, analogous to the standard odd-even mass 
staggering caused by T = 1 pairing. Although the T=0 ground state in odd-odd nuclei is of 2qp character, it is still 
an open problem of how to construct a proper trail wave function appropriate for odd-odd nuclei. It seems that the 
isospin restoration (see discussion in Ref. p^) may play a key role here. 



C. The cranking calculations 

As mentioned previously, our schematic, separable np-interaction does not allow for T=0 act np-pairing. In a 
condensate, that is dominated by T=0 aa pairs, there is little resistance for these pairs to align their angular 
momenta. Note that this alignment occurs without invoking any pair-breaking mechanism. The alignment is a 
smooth function of the frequency, and no backbend occurs, see Figure 3 in our letter [^ . 

In the low spin regime, this is of course an unphysical behaviour. The presence of T=0 Ckct pairing will invoke 
the pair-breaking mechanism (low- J np-pairs dominate at low-spin [ll|] ). However, as long as we do not use 
any physical measure sensitive to the isospin of the np-pairs the T=0 aa np-pairing can be mocked up by T=l aa 
np-pairing and vice- versa via a simple readjustment of the strength (isospin-broken model). Indeed, both modes will 
respond to nuclear rotation in similar way, i.e. through the pair breaking mechanism which leads to a suppression of 
this field at high spins. 

In contrast to the aa pairing the aa type pairing will become enhanced at high rotational frequencies as Coriolis 
and centrifugal forces align more and more pairs along the axis of collective rotation. Indeed, as shown in |p^ , po| , the 
T=0 phase composed of high- J np-pairs survives to much higher angular momentum. This is in accordance to recent 
shell model Monte-Carlo calculations . At very high spins, our schematic T=0 aa pairing is probably somewhat 
too strong since there is no restriction in the J-values of the pairs that participate in the scattering. However, in 
accordance to early calculations, we expect to pick up the main essence of the high spin behaviour. 

Fig. shows an illustrative example phase-diagram of the critical frequency (the frequency corresponding to the 
onset of np-pairing of aa) versus the strength of this force. The calculations were performed using the BCSLN model 
for the N — Z nucleus ^^Cr at constant deformation and for huj < 2MeV. Filled circles represent calculations with 
Gpp = Gnn = G"p while open circles correspond to an isospin symmetry broken Hamiltonian with G"p = 1.3Gpp(^nny 
Two things are worth to be noticed: (i) Tiujcrit is very sensitive to the strength parameter and (ii) the alignment (in 
this case the alignment of fi /2 quasiparticles) triggers the onset of T=0 np-pairing setting a natural threshold (x"") 
for the process. Above this threshold i.e. for a;"" < x"" the value of hujcrit increases very sharply. Therefore, even 
the optimistic scenario would involve the onset of rtp-pairing (if at all possible) around the crossing frequency i.e. in 
the region which is anyhow the most difficult to describe theoretically. This is in accordance to single-j shell model 
calculations 
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D. Terminating states 



In the p/-shell, the state of the art Ohuj sheU-models are able to describe nuclear structure with excellent accuracy 
up to the maximum spin states (terminating states) which is built upon the pure /7/2 sub-shell, see e.g. [|9|. For 
example, for the case of ^*Cr, the state of maximum spin of / = W^Ti can be reached within irfy^^i/fy^ configuration. 
Particle-hole excitations to and sub- shells allow further to built states of spin up to / = 20+ ?i within the 
full pf space. However, already at these spins the states involving cross-shell particle-hole excitations are expected to 
compete. Unfortunately, with the present-day shell-model techniques it is not feasible to perform calculations within 
the model space needed to describe states above spin / — 16h, which in principle should include the whole srf-shell, 
p/-shell and the gg/2 sub-shell. 

In the standard cranked mean- field model calculations the collectivity of the rotational band in ^^Cr is essentially 
exhausted after the simultaneous alignment of fi/2 neutrons and protons p9| and one enters the non-collective, 
unpaired regime around Ix = 16h. Building up higher angular momentum states proceed further in discrete jumps, 
which are due to crossings between single-particle down-sloping and up-sloping routhians at spherical shape. This 
reoccupation process will finally lead again to the onset of collectivity. 

The presence of T=0 correlations at high angular momenta may strongly alter the behaviour of high spin states 
discussed above. In our model, above the terminating / — 16h state we observe smooth increase of Ix versus huj 
rather than the step-like process expected in standard scenario, see Fig. ||. This is due to the pair scattering from 
^3/2 and /7/2 into the aligned .gg/2 and f^^2 orbits which is entirely due to T=0 pairing. Partial occupation of these 
orbits triggers the onset of collectivity after the terminating state [^. New experiments targeted to measure the 
evolution of the rotational band in ''^Cr beyond / = 16?i may therefore shed light into the nature of the mean-field 
T=0 pairing correlations at high spins. Particularly, the possibility of the T=0 pairing correlations to affect the 
electromagnetic transition rates. Indeed, without T=0 correlations, the wave functions of high-spin states will have 
the dominant component of 2p-2h and/or 4p-4h configuration (due to promotion of a 7^/7/21^/7/2 pa-ii' or 1^/7/2^/7/2 
quartet to higher sub-shells). The E2-decay from such a configuration into the Op-Oh i'n'/y2 ® ^17/2) state would then 
be strongly hindered. In contrast, T=0 pairing will cause mixing of these (multi) particle- (multi) hole configurations 
(collectivity) enhancing the E2 transition amplitudes. 



E. TRS-calculations at superdeformed shape 

Nuclear superdeformation (SD) is an extreme case of the spontaneous breaking of spherical symmetry. In this case 
our T=0 pairing interaction, scattering uniformly np-pairs, may even be considered as quite a realistic approximation 
provided that the strength is properly adjusted. The standard total routhian surface (TRS) as well as Skyrme- 
Hartrec-Fock calculations predict SD bands to appear at relatively low excitation energy in TV ~ Z nuclei of A«80-90. 
Predicted deformations are /32 ~ 0.55 for ^^Ru and even as large as /32 ~ 0.75 for *^Mo, see [ |50| , |5]| . To illustrate 
the influence of T=0 interaction on nuclear superdeformation we performed TRS calculations for SD ^^Ru. Three 
variants of the calculations different in the treatment of pairing have been performed: (i) unpaired single-particle {ii) 
including isovector T=l pairing only, {Hi) including both T=0 and T=l pairing with x^^'^ = 0.9, x'^^^ = 1.1, and 
x^^'^ = 1.3. The dynamical moments of inertia (J^^-* = dix/duj) versus rotational frequency resulting from this set of 
calculation are shown in Fig. 

As expected, the calculations (i) and (ii) differ strongly at low spins but converge at higher frequencies where static 
T=l pairing correlations are essentially washed out by the Coriolis force. The pronounced kink at huj « 0.75 MeV 
in calculations (ii) is due to the simultaneous alignment of the /iii/2 protons and neutrons. Choosing the strongly 
undercritical strength of G^^^ , corresponding to = 0.9, results of course in a pair field dominated by T=l pairs. 
However, after the /iii/2 alignment and the disappearance of static T=l pairing correlations, the T=0 correlations 
start to build up, resulting in a somewhat larger moment of inertia as compared to the unpaired case. Clearly, the 
onset of T=0 pairing is triggered by the alignment. The x'^^^ = 1.1 corresponds to a slightly undercritical G"^^" 
strength at low huj. The low frequency hump in J^^^ at Tiio « 0.3 MeV is in this case due to the rotation-induced 
phase transition from predominantly T=l to T=0 pairing p8|j2(|_. For larger values of the hump moves towards 



lower frequencies in accordance with the phase diagram in Fig7]7| discussed in Subsect. VC, and the /ij^^/o alignment 



is smoothed out. Moreover, the high-frequency part of J^^^ increases exceeding the 'rigid body' value □ by up to 



We use the notation of 'rigid body' value as the one given by calculations without pairing. Note that this moment of inertia 
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Ri30%. Although the increase depends very much on the strength of the T=0 interaction (and cssentiaUy vanishes for 
x^^°<0.8) the observation of systematic deviations from the standard T=l model predictions can be considered as a 
rather robust indicator of the importance of the stretched i.e. oca type T=0 correlations at high spins. 

VI. SUMMARY AND CONCLUSIONS 

Our work presents a number-projected, cranked mean-field formalism appropriate for the simultaneous treatment 
of T=l and T=0 correlations. The number-projection technique is a generalization of the Lipkin-Nogami method. 
We use a schematic np-pairing interaction because our interest is to study general properties of np-paired solutions 
rather than to reproduce properties of specific nuclei. 

Our work shows that the cxclusivcncss of different pairing phases is inherent to the BCS method. Already ap- 
proximate number-projection leads to a mixing of the T=0 and T=l pairing phases. Common for both approaches, 
though, is that different pairing phases counteract each other. 

It is also demonstrated that neutron or proton excess block np-pairing. This generalized blocking mechanism causes 
np-pairing to vanish rapidly when departing from the N = Z line. The calculations suggest that already at \Tz \ ~ 2 
there is essentially no chance to observe any collective np-pairing. This mechanism serves as a possible microscopic 
explanation of the Wigner energy at the level of mean-field theory. The even-even versus odd-odd mass staggering 
in T=0 states oi N — Z nuclei indicates a 2qp structure (one broken np-pair) of the T=0 states in odd-odd N = Z 
nuclei. This will require blocked calculations which are beyond the scope of our present paper. 

Our model predicts that the T=0 correlations do not vanish at high spin, in accordance to previous calculations. 
In the high-spins regime stretched T=0 pairing is predicted to be dominant. The onset of this pairing phase seems 
to be triggered by the alignment of high-j quasiparticlcs and the related quenching of isovector pairing correlations. 
The presence of T=0 pairing at high spins is predicted to increase the moments of inertia as well as to affect the 
evolution of rotational bands beyond their standard terminating states. The size of this effect depends critically on 
the strength of the T=0 force in the aa channel. Below a certain limit, in our case below x'^^<0.8, the effect will 
not be seen. 

This research was supported in part by the U.S. Department of Energy under Contract Nos. DE-FG02-96ER40963 
(University of Tennessee) , DE-FG05-87ER40361 (Joint Institute for Heavy Ion Research), DE-AC05-96OR22464 with 
Lockheed Martin Energy Research Corp. (Oak Ridge National Laboratory), by the Polish Committee for Scientific 
Research (KBN) under Contract No. 2 P03B 040 14, and by the Swedish Institute. 



strongly differs from that of a rigid body, which would be a constant. 
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APPENDIX A: THE TIME-REVERSAL INVARIANCE, ISOSPACE DEFORMATIONS AND THE BCS 

SOLUTIONS. 



The time-reversal invariance implies the following relations for the hermitian density matrix and antisymmetric 
pairing tensor elements |T^]: 



< 



PaT,aT = Par,aT 
PaT,aT' — P-^r,7iiT' 
Par^'ar 
^ PaT,a-T = ^PaT,a-T , 



K. 



a.r.CL — T — ' 



* 



(Al) 



The standard BCS constraints (AiV) = (AZ) = automatically imply that (T^) — {N — Z)/2. Superimposing 



further symmetry constraint (T.^) — (Tj,) = together with time-reversal symmetry induced relations (Al) implies: 

(Ta;) = ^ Re(pap_a„) = (T^^) = ^ Im(pap,Qn) = /9ar,Q-T=0 V u. (A2) 

a>0 a>0 

The HFB solutions further demand that the generalized density matrix: 

p K 



3ft 



K' 1 p 



(A3) 



obeys the idempotency condition 5?^ — 3? or, equivalently, satisfies the following relations: pn = up* and p^ p = 
KK* . These relations impose six constraints on the density matrix and pairing tensor elements which all must be 
fulfilled simultaneously. They can be written as: 



^aT,a—T PaT,a — T ^" ^ar ,a — T Par .a — r 



(A4) 



r,OL — rh'o/_r,cx—T} 



(A5) 



(A6) 



^aT^a—ri^PaT^ar Pa — T-a — r') PaT,a — T(^^a.T.Oi.T ^a — T,Q. — T^ 



Par ar ~^ |/^ar,Q — r| ~h |par,a— r| Par^ar 



^aT,a. — T 



(A7) 
(A8) 



Par,a— r(l Pcxt^olt Pa — T,OL — T^ ^a.T,a — Ti,^aT,a.T ~l~ '^a — r,Q — r) 



(A9) 



Parrot— T ( 1 Par.ar Pa — r.a — r^ 



(AlO) 



For N ~ Z nuclei, the coherent np-paired solution can be obtained when Pap,ap — Pan,an and n^pj^p = —Kan.an- 
This is the solution found by Goodman and we refer reader to the Ref. |Q for f urth er det ails. For N ^ Z nucleus 
however, neither |pQp,Qp| = |pan,an| nor |Kctp,Qp| = Using the relations (A4)-(A1C) it is straightforward to 

show that the additional symmetry (axiality in isospace) ( [A^ ) rules out np-pairing of olol type. Indeed, Eqs. (A7) 
and (A9) give instantly k^tR-t = 0- 
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FIG. 1. Value of the a) Lipkin-Nogami parameters A'^^ and b) contributions to the total energy SE^^ff), calculated for the 
= nucleus **Cr and = 1 nucleus ^"Cr. The total correction to the binding energy (filled triangles) and the constituents 
due to proton/neutron pairing (O) and np-pairing (o) are shown separately. Note, that the contribution due to np-pairing is 
effectively repulsive. 

FIG. 2. Schematic representation of the generalized blocking mechanism due to the neutron excess. Shaded area indicates 
the single-particle levels which do not contribute to np-pair scattering. 

FIG. 3. Experimental (•) and calculated strength W^(/l) of the Wigncr energy for p/-shcll nuclei. The curves marked as (o) 
correspond to the BCSLN model calculations with x^^" = 1 and 1.30, respectively. The curves labeled by (A) denote the BCS 
model calculations with a;^^" = 1, 1.13 and 1.15, respectively. 

FIG. 4. Values of x'^="/xj=l^ versus T^. Open circles preseut the BCS model results {xq " « 1.13) while solid circles mark 
the BCSLN calculations (xq^^ « 1.30) in Cr nuclei. The line denoted by (★) represents BCSLN calculations in Sr nuclei 
(xq^" « 1.25). See text for further details. 

FIG. 5. Calculated mass excess AE = E{x^^°) — E{x^^'' = 1) as a function of Tz ■ The calculations have been performed for 
Cr (open symbols) and Sr (filled symbols) isotopes and for a;^^° = 1.1 (A), 1.2 (O), 1.3(o). The line denoted by (★) corresponds 
to the 'physical' value of a;^^" = 1.25 representative for Sr isotopes. The strength parameters for like-particle correlations were 
calculated using the average gap method. In our case it yields Gpp « 0.385(0.250) MeV for Z=24(38), respectively. The 
difference in the pairing strength, AG = G„„ — Gpp versus is shown in the insert. 



FIG. 6. Experimental (•) and calculated (o. A) correction to the Wigner energy, d{A)T=o- Note, that both BCS and BCSLN 
theories give essentially no contribution. For further detail see text. 

FIG. 7. Frequency Tuocru for the onset of the T=0 aa np-pairing as a function of the strength a;^"" of this force. Calculations 
were performed for **Cr at constant deformation (32 = 0.25 using the BCSLN model. The line marked by (•) shows the 
calculations with an isospin-symmetric hamiltonian i.e. for Gj^^ = Gnn^ = G^p^ while (o) denotes the calculations for an 
isospin- broken hamiltonian with G^j^^ = 1.3Gpp(„,i). The vertical arrow marks the position of the crossing frequency of the 
fr/2 quasiparticles calculated without np-pairing. 

FIG. 8. Gain in angular momentum, Six, due to the onset of T=0 aa np-pairing relative to the predictions of a model 
without np-pairing. The quantity Six was calculated for **Cr {^2 = 0.25) at high frequency, fiu) = 2 MeV and is shown as a 
function of the strength parameter a:^^°. 



FIG. 9. The dynamical moment of inertia for the SD band in the N = Z rmclcus ***Ru as a function of huj. Unpaired (T=l 
paired) calculations are marked with o (•), respectively. The calculations with T=0 pairing are labeled by A (a;'^~°=0.9), □ 
(a;^^°=l.l), and O (a;^^°=1.3). Note the large moments of inertia for the T=0 paired solutions at high frequencies. 
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